Local contact numbers in two dimensional packings of frictional disks 
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We analyze the local structure of two dimensional packings of frictional disks numerically. We 
focus on the fractions Xi of particles that are in contact with i neighbors, and systematically vary 
the confining pressure p and friction coefficient p,. We find that for all p,, the fractions Xi exhibit 
powerlaw scaling with p, which allows us to obtain an accurate estimate for Xi at zero pressure. 
We uncover how these zero pressure fractions Xi vary with (i, and introduce a simple model that 
captures most of this variation. We also probe the correlations between the contact numbers of 
neighboring particles. 

PACS numbers: 45.70.-n, 46.65.+g, 83.80.Fg 



While soft frictionless spheres experience a critical jam- 
ming transition in the limit of zero pressure, where prop- 
erties such as elastic moduli, contact number, density, 
characteristic frequencies and lengthscales exhibit pow- 
erlaw scaling P, H H, 0| j the situation is more delicate for 
frictional systems. The approach to the jamming tran- 
sition is still governed by the pressure, p, but a range 
of densities and packing properties can exist depending 
on the value of the friction coefficient p, the mobiliza- 
tion (ratio of frictional to normal forces) of the frictional 
contacts and the packing history [E, S 0, S| • In partic- 
ular, in d dimensions, the contact number at jamming, 
z c , can take on a range of values between d + 1 and 2d, 
in contrast to frictionless sphere packings which always 
reach their respective isostatic contact number z? so = 2c? 
at jamming. The proximity to the isostatic contact num- 
ber governs the scaling near jamming — for frictionless 
spheres, properties such as elastic moduli scale with dis- 
tance to jamming. However, for frictional packings these 
properties only scale with distance to the isostatic limit 
Zjl == d+ 1, and in general not with distance to jamming 
[if M 0] j although this depends on whether fully mobi- 
lized contacts are treated as frictional or slipping [§]. 

We recently studied the case of frictional spherical 
disks in two dimensions, and focussed on packings that 
were equilibrated very gently @, 0, Q . This eliminates 
preparation history and mobilization as unknowns: for 
given pressure p and friction coefficient //, packings with 
well defined statistics are obtained. The gentle equilibra- 
tion procedure also allows to approach the isostatic limit 
for frictional systems, z c = z^ so = d+ 1 when p, — > oo and 
p —> — here jamming has many of the critical features 
observed for frictionless systems [1, 0] • 

One additional surprise is that for finite values of /x, 
such gently equilibrated packings still reach a generalized 
isostatic limit @, @|. This means, in short, that a sub- 
stantial number of contacts get fully mobilized, i.e., their 
frictional forces ft satisfy the bound |/ ( | < /„, where /„ 
denotes the normal force. If these fully mobilized con- 
tacts are seen as slipping, the critical nature of the vi- 
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FIG. 1: Variation of the fractions Xi(p = 0, p) of particles with 
i — 2, 3, . . . , 6 contact neighbors as function of the friction 
coefficient p. The full curves are predictions from a simple 
model (Eqs. [T]) with fixed variance a 2 = 0.6. 



brational density of states at jamming is restored for all 
values of \i [§]. 

Here we probe the fractions Xi (p) of particles that have 
i contacts for these frictional packings. These fractions 
are the simplest characteristics of the contact network 
beyond the average contact number z. It is thus natural 
to ask how the fractions x% depend on p and /i. We find 
that, for given /i, the fractions Xi(p) exhibit scaling with 
p similar to the scaling of the total contact number z. 
This allows us to extrapolate these fractions to p — > 0, 
and this is the case on which we focus our attention. As 
is shown in Fig. [TJ the fractions Xi vary substantially with 
/i, and reach well-defined values in the limits where p — > 
or /i — > oo. We find a number of simple but unexpected 
relations between the various Xi, and introduce a simple 
model that, given z(n), gives a good prediction for Xi(/i). 

Packings — Following the numerical systems un- 
der consideration are two dimensional packings of 1000 
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FIG. 2: Contact fractions function of pressure, (a-c) 

For three representative values of /i, the Xi scale linearly with 
p 1 ^ 3 (equivalent to </> 1/ ' 2 for the hertzian interaction), and we 
are able to extrapolate to p — 0. (d) The sum £3 + x$ ~ 0.5, 
for all values of jj, and p studied. 



spheres with 20% polydispersity in the diameter of the 
particles in a square box with periodic boundary con- 
ditions. The grains interact through 3c? Hertz-Mindlin 

forces, i.e with the normal force /y between particles i 
3/2 

and j proportional to 8^ , with 6ij the overlap of the 
two particles. The Young modulus of the grains is set to 
1, which determines the pressure unit, and the Poisson 
ratio is set to zero, while the unit of length is the aver- 
age grain diameter. The construction and equilibration of 
the packings has been described in detail elsewhere 0, Q • 
Rattlers, particles which have no appreciable interactions 
with any of the other particles, are always left out of the 
analysis of the packings and contact statistics. For each 
value of /1 e [1(T 3 , 10 3 ] and p e (lfr 6 , 1CT 3 ), 30 config- 
urations were generated independently. 

Scaling of fractions Xi with pressure — As is shown 
in Fig. [2^-c, Xi(p, /1) scales linearly with p 1 ^ 3 , which al- 
lows us to extrapolate their values for finite p to the 
(un)jamming limit at p = 0. This scaling is the same 
as the scaling of the total contact number z with p, 
which for the Hertzian interactions employed here is con- 
sistent with the scaling that the excess contact number 
Az := z — z c scales with the square-root of the excess 
packing fraction. This relation is well known for friction- 
less systems Ql , but also appears to hold for frictional 
systems [1, [Tl[ — our data here suggests that it also 
holds for the individual contact fractions, irrespective of 
the value of p.. 

A second robust finding is illustrated in Fig. [2ji: the 
number of particles that have an odd number of con- 
12j, is close to 1/2 — the number of particles with 



tacts 

an even or odd number of contacts is therefore approxi- 



FIG. 3: Rate equation models for the equilibrium contact 
fractions. (a,b) A model with a single rate u> is sufficient for 
/1 — > and y, — * 0. (c) For finite u, we introduce individual 
rates un which correspond to the total decay rate for contact 
number i. 



mately equal, irrespective of pressure or value of [i. We 
do not have a satisfactory explanation for this. 

The extrapolated fractions Xi at jamming — In the re- 
mainder of this paper we focus on #*(//) at zero pressure. 
Since Xi has to be zero for i = 1, and the fraction of 
particles with 7 contact are negligible for the polydisper- 
sities employed here we focus on i ranging from 2 to 6. 
As shown in Fig. [TJ the variation of Xi with \i is greatest 
for /1 between 0.1 and 1, with the small and large /i limits 
apparently well behaved. 

The functional forms of Xi(fj.) for i—3 and 5 are similar, 
as are the functional forms of for i—2 and 4. This is 
related to the observations that xz + X5 « 1/2. One also 
notices that, approximately, x n (fi — > 0) ss x n +i{^ — > 00). 
In fact, for small /Lt, the fractions x% and x§ tend to 1/4, 
while Xi approaches 1/2 — for large /z, X2 and X4 tend 
to 1/4, while £3 approaches 1/2. 

In the limits /x = or fj, = 00, we can estimate these 
fractions by a very simple argument. Let us first focus on 
the zero friction case. Assuming that there are only par- 
ticles with three, four or five contacts, the fractions £3, X4 
and X5 can immediately be calculated, since combining 
the condition that X3 + X4 + x$ — 1 with the isostaticity 
condition 3a;3 + 4x4 + 5x§ — 4 implies X3 — x§ , and hence 
X3 = 1/4, .T4 = 1/2 and x$ — 1/4 — a similar argument 
holds for X2,X3 and X4 in the limit of infinite friction. 
Deviations from this result arise since a small fraction 
of particles with respectively six and five contacts arise, 
weakly breaking the "three particle species " condition 
underlying this argument (see Fig. 1). 

Simple rate equation model — The ratios X3/X4 = 
X5/X4 = 1/2 can also be understood in terms of a simple 
stochastic model where we imagine distorting a certain 
packing, creating and breaking contacts but keeping the 
overall contact number and the ratios Xi constant. In 
the case of three species only, particles with 4 contacts 
can become 3's and 5's, while 3's and 5's can only be- 
come 4's (See Fig. 3). Since the transition probabilities 
must all be equal (since always two particles take place in 
such an event), and, on average, we require the fractions 
Xi to be constant, we get, in this simple approximation, 
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= 2^3 = 2x5. This heuristic argument can be written 
as a rate equation model, as shown in Figure^. Once we 
normalize the rates such that the total decay rate of each 
species is 2u>, we obtain as steady state £4 = 2x 3 = 2x 5 . 

For intermediate values of /i, the number of species is 
four (if we neglect a small number of z = 6-contacts). A 
single decay rate would than imply that {x 2 ~ 1/6, £3 w 
1/3, £4 w 1/3, £5 « 1/6} and z — 3.5 — clearly a single 
rate does not capture the data. Figure [3fc shows an ex- 
tended model where we now associate an individual rate 
u>i to each species i, so that the total decay rate of that 
species is 2loi. The solution to this model is Xi ~ 1/uJi 
for i = 3, 4 and Xi ~ 1/ (2u>j) for z = 2, 5. 

Explicit solutions of rate equation model — We now 
seek an explicit solution of the four species model for 
the contact fractions as a function of the friction coeffi- 
cient. To achieve this, we introduce two constraints on 
the model beyond the trivial normalization constraints 
Si=2 Xi ~ 1 an d Et2 = z {lA- Fhst, we constrain 
our model by the empirical observation that the num- 
ber of particles with odd and even contacts is equal, i.e., 
£3 + ^5 = 0-5- Additionally, we impose the variance of 
the contact fraction distribution, ^2 i=2 x i{ z ~ *) 2 = 0-2 • 
The solution to the resulting set of equations is 

x 2 = ({z - 4) 2 + a 2 - 1/2) /4 

x 3 = (-(z-3) 2 -a 2 + 5/2)/4 (1) 

x A = (-(z-4) 2 -a 2 + 5/2) /4 

x 5 = ({z - 3) 2 + a 2 - 1/2) /4 

To obtain definite predictions from this set of equa- 
tions, we need to determine the variance a 2 . In the ex- 
treme limits, and under the simplifying assumption that 
only three species with fractions 1/4,1/2,1/4 arise, we 
find a 2 — 0.5 (notice if more species are present, a 2 will 
be larger). Fixing now a 2 = 0.5 over the whole range 
of friction coefficient, we obtain the prediction shown in 
figure |4] There are no additional fit parameters to this 
solution, and the agreement is quite good. 

We have numerically studied the actual variance of a 2 
from the data, and find that for our data it varies between 
0.57 and 0.65 — when we fix a 2 — 0.6, the fit becomes 
significantly improved, as shown in Fig. 1. 

Correlations — The rate equation model derives from 
its implicit assumption that the contact numbers of par- 
ticles and their neighbors are uncorrelated. Based on this 
assumption, we can calculate the theoretical fraction q\^ 
of contacts between particles with i and j, given Xi and 
Xj . Since the total fraction of contacts for particles with 
i contacts is given by ixi/z, the uncorrelated prediction 
for qij is 

^ = 2w iforMi . $ = Wp_ fa[ i = j (2) 

Figured shows the ratio qij/qlj of the observed frac- 
tion of contacts and the uncorrelated prediction [l3j . For 




FIG. 4: Contact fractions as a function of /i in the extrapo- 
lated limit p — > 0. The curves show the model solution from 
equation[TJ with a variance a 2 — 0.5. 

intermediate values of i and j the prediction is quite rea- 
sonable, as qij/qfj remains bounded between 0.5 and 1.6 
orso. Contact pairs with very dissimilar i and j are fa- 
vored - this is likely an effect of polydispersity, since small 
particles with few contacts prefer to sit next to larger par- 
ticles with more contacts. A detailed study of this is left 
for the future. 

In Figure [5]b we have divided out the ratio at an inter- 
mediate /i, to more clearly see the variation of q^ with /z. 
This shows that the fractions corresponding to particles 
with Xi that are abundant (such as 944 for small fj, and (733 
for large /i) do not vary strongly with fi. There appears 
to be a correlation between the relative over representa- 
tion of contacts and the overabundance of the species of 
particles (i.e., for large \x, there are many particles with 
2 or 3 contacts, and (723 is over abundant, while there are 
very few particles with 4 and 5 contacts, and the ratios 
944, 945 and 955 are even less likely) — we have no clear 
explanation for this. 

Outlook — Simple arguments allow us to estimate the 
contact fractions Xi, which can be seen as fingerprints of 
the system. Since frictional systems depend on history, 
we expect the fractions and their variation to be a useful 
step in identifying the effects of preparation history be- 
yond average values such as overall contact number and 
density. 
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